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We use a Monte Carlo approach to study hadron azimuthal angular correlations in high energy 
proton-proton and central nucleus-nucleus collisions at the BNL Relativistic Heavy Ion Collider 
(RHIC) energies at mid-rapidity. We build a hadron event generator that incorporates the produc- 
tion of 2 — > 2 and 2 — > 3 parton processes and their evolution into hadron states. For nucleus-nucleus 
collisions we include the effect of parton energy loss in the Quark-Gluon Plasma using a modified 
fragmentation function approach. In the presence of the medium, for the case when three partons 
are produced in the hard scattering, we analyze the Monte Carlo sample in parton and hadron 
momentum bins to reconstruct the angular correlations. We characterize this sample by the number 
of partons that are able to hadronize by fragmentation within the selected bins. In the nuclear 
environment the model allows hadronization by fragmentation only for partons with momentum 
above a threshold py loah = 2.4 GeV. We argue that one should treat properly the effect of those 
partons with momentum below the threshold, since their interaction with the medium may lead 
to showers of low momentum hadrons along the direction of motion of the original partons as the 
medium becomes diluted. 

PACS numbers: 25.75.-q, 25.75.Gz, 12.38.Bx 
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I. INTRODUCTION 



One of the main experimental discoveries in the field of 
high energy heavy-ion reactions has been the suppression 
of particles, with momentum similar to the leadingone, 
in the away side of azimuthal angular correlations |l| . It 
is by now believed that this phenomenon is caused by 
the energy loss of partons moving through the strongly 
interacting medium formed in the aftermath of the re- 
action. This picture is far from being that simple since 
when the difference in momentum between leading and 
away side particles increases, either a double peak struc- 
ture or a broadening of the away-side peak appears, while 
cither of these are absent in p + p collisions at the same 
energies Q. 

Given that the structures in the away side are more 
prominent for relatively small momentum particles, the 
above features have generated different explanations 
based on medium effects. The current trend explains the 
double peak/broadening in the away side as due to initial 
state fluctuations of the matter density in the colliding 
nuclei. These fluctuations would in turn be responsible 
for an anisotropic flow of partially equilibrated low mo- 
mentum particles with the bulk medium. Recent descrip- 
tions of this scenario, based on hydrodynamics, have re- 
produced successfully the experimental V3 Q. Nonethe- 
less, it has also been shown recently that there is a strong 
connection between the medium's path length and the 
observed away side structures [|[ . This connection is ex- 



pressed through the dependence of the azimuthal corre- 
lation on the trigger particle direction with respect to 
the event plane in such a way that, for selected trig- 
ger and associated particle momenta, the double peak is 
present /absent for out of plane/in plane trigger particle 
direction. 

The difference in path lengths traveled by partons in 
the medium is at the core of the idea that the away side 
structures may include contributions of 2 —> 3 processes. 
For instance, when the hard scattering resulting in this 
kind of processes happens close to the edge, there is a 
large chance that at least one of the three final state par- 
tons travels a large distance through the medium there- 
fore, having a high probability of being absorbed. The 
resulting hadronic process has two particles in the final 
state but its origin is an underlying partonic event with 
three partons in the final state, one of which was unable 
to hadronize. Conservation of momentum at the parton 
level followed by collinear fragmentation gives rise to a 
distinctive angular dependence in the azimuthal correla- 
tion whereby, the angular difference between leading and 
away side particles is close to 2ir/3 radians, regardless of 
whether one or two partons in the away side survive their 
passing through the medium and are able to hadronize. 

The above scenario has been put forward and explored 
in Refs. [f|, using the leading order QCD matrix elements 
for 2 — > 3 parton processes. A serious limitation of this 
approach stems from the way the collinear singularities 
are avoided. This was implemented by restricting the 
phase space for parton production to the regions where 
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the angular difference between leading and away side par- 
tons is far from and ir. To overcome such limitation, in 
this work we present a Monte Carlo approach to study 
azimuthal correlations in p + p and A + A collisions. We 
use the MadGraph5 Q event generator, which includes 
built-in functions to cancel collinear and soft divergences 
for 2 — > 3 parton processes. To study medium effects we 
use the modified fragmentation function approach. We 
also study 2 — > 3 parton processes that contribute to only 
two hadrons in the final state. The work is organized as 
follows: In Sec.|TT]we present the basics of the description 
for three hadron production in p+p and A + A collisions. 
In Sec. IIIII we introduce a Monte Carlo event generator 
to implement the calculation of azimuthal angular corre- 
lations. To study the details of the effect of energy loss 
on partons within different momentum bins, in Sec. IIVI 
we present both the parton and hadron px distributions 
for the Monte Carlo generated events. We note that for 
the model parameters used, low px hadrons come from 
relatively high pt partons. We use these samples to build 
the azimuthal angular correlations in Sec. El Finally, we 
discuss our results and conclude in Sec. IVI1 in particular 
we argue that the double peak/broadening in the away 
side observed in data may still be described within this 
approach, by considering the interaction of those partons 
that are not able to hadronize outside the medium, but 
still interact with the bulk partons. 

II. THREE-HADRON PRODUCTION 

In previous works [1[ we have computed the differen- 
tial cross section for three-hadron production in p + p 



and A + A collisions at mid-rapidity. In order to show 
the essentials of that approach, we hereby summarize it 
and refer the reader to the aforementioned references for 
further details. 



In the case of p + p collisions the differential cross 
section was obtained by convoluting over the incoming 
momentum fraction, the initial distributions / p for par- 
tons within the colliding protons, the matrix elements 
|^W 2 ^ 3 1 2 and the fragmentation functions D F / H . For this 
we used the CTEQ6 parametrization Q, the leading or- 
der matrix elements describing the process at the parton 
level [|[ and KKP fragmentation functions , for a given 
total center of mass energy yfs available for the collision. 
We considered collinear fragmentation thus, the angles 
that define the direction of the away-side hadrons, Of 
(i = 2,3), are linearly related to the the parton angles Oj 
= 2,3). 



To consider the process within a central heavy-ion col- 
lision and thus account for the effects of energy loss, we 
resorted to the model put forward in Ref. ■ The model 
considers an initial gluon density obtained from the over- 
lap of two colliding nuclei, each with a Woods-Saxon den- 
sity profile. The gluon density of the medium is diluted 
only due to longitudinal expansion of the plasma since 
transverse expansion is neglected. 

The gluon density p g is related to the nuclear geometry 
of the produced medium. We use the modified fragmen- 
tation functions M 
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where z' nm = h m / (p n — AE n ) is the rescaled momentum 
fraction, of hadron H m originated from the fragmenting 
parton P„, z' nm . g = (L/ X)(h m / AE n ) is the rescaled mo- 
mentum fraction of the radiated gluon, AE n is the aver- 
age radiative parton energy loss and (L / X) is the average 
number of scatterings. Both AE n and (L/X) are related 
to the gluon density of the produced medium. This re- 
lation comes in by considering the path length traveled 
through the medium by the partons produced in the hard 
scattering. The path length is determined from the ra- 
dial position r and the angle </> that the leading parton 
makes with the radial direction, within the medium. 

The average energy loss AE n is also proportional to the 
one dimensional energy loss {dE n / dL)id parametrized as 





- 1.6 


1.2 


7.5 H 


_ Pit 


-l 




-/J 











The one-dimensional energy loss per unit length parame- 
ter, e , is related to the mean free path Ao by eoAo = 0.5 
GeV. Notice that Eq. (2) means that there is a lower 
threshold for p n given by p!£ roah = 1.6 x p = 2.4 GeV 
below which the one-dimensional energy loss becomes 
meaningless. This is interpreted as reflecting the fact 
that the energy loss model cannot produce hadrons by 
fragmentation for partons below this momentum thresh- 
old. Those partons should thus be absorbed by the 
medium and their effect has to be accounted for by other 
hadronization mechanisms. 

The above described procedure has an inherent limita- 
tion which we have already identified in Refs. Q. These 
are related to the angular cuts required to avoid the 
collinear singularities when two of the final state hadrons 
are either parallel or antiparallel. In order to overcome 



3 



such limitations we proceed to present a Monte Carlo 
implementation of two and three hadron production in 
p + p and A + A environments dedicated to the study of 
azimuthal angular correlations. 



III. MONTE CARLO EVENT GENERATOR 
WITH ENERGY LOSS 

A. Parton level event generator 

We implement the production of parton events using 
MadGraph5 @ , in p + p collisions at RHIC energies at 
mid-rapidity (|?7| < 0.5). MadGraph5 generates the he- 
licity amplitude subroutines to build the hard scattering 
matrix elements, used as inputs for the partonic event 
generator. It also generates a set of functions, inspired 
by the dipolc subtraction method [ll|, that remove the 
soft and collinear singularities in the phase-space inte- 
grated matrix elements arising under certain kinematical 
conditions, producing an unweighted set of hard scatter- 
ing partonic events. A typical run involves declaring the 
values for certain parameters such as the beam energy, 
the factorization scale /i 2 , the fragmentation scale Q 2 , 
the px bin size of the final state partons and the value of 
parameters like xqcut that controls the numerical can- 
cellation of the singularities, according to the particular 
kinematics. It is worth mentioning that, as suggested in 
the literature [f|, we use values for xqcut below 10 GeV. 

B. Hadron level event generator 

In order to achieve hadron events we take the parton 
level events and evolve the partons into hadrons with 
collinear fragmentation. 

To study p + p collisions, partons are fragmented into 
hadrons, by means of KKP fragmentation Q. To de- 
scribe medium effects on the propagation of partons and 
their hadronization in a nuclear collision, we include en- 
ergy loss effects through KKP modified fragmentation 
functions fioj . In the former case, parton i with momen- 
tum p is assigned a random value of momentum fraction 
z, which allows it to evolve into hadron H with momen- 
tum p H . In the latter, we take into account that the hard 
scattering can take place anywhere within the nuclear 
collision overlap area and that the produced partons can 
travel in any direction. This is implemented by assigning 
the hard scattering a random radial position r and a ran- 
dom direction (j) that the leading parton makes with the 
outward radial direction. These parameters are used to 
determine the average number of scatterings (L/A), the 
energy loss AEi for each parton and their modified mo- 
mentum fractions z\. In this way the hadron level events 
are characterized by the particle species and their kine- 
matical variables both in p+p and A+A collisions. Notice 
that the KKP fragmentation functions are valid only for 
hadron momentum fractions z in the range 0.1 < z < 0.8 



which affects the shape of the azimuthal correlations for 
the away-side hadrons' relative angular separation close 
to and it radians. 

For the nuclear environment, since the minimum al- 
lowed parton momentum for hadronization is, according 
to the model in Ref. [h|, px = 2.4 GeV, then we have 
three different types of events corresponding to the num- 
ber of partons that arc absorbed by the medium due to 
energy loss. Notice that this implies that there is a con- 
tribution from 2 — > 3 partonic processes (where one of 
the final state parton is absorbed by the medium) to the 
di-hadron azimuthal correlation function. 



IV. p T DISTRIBUTIONS FOR 2 -> {2, 3} 
HADRON PRODUCTION 

In order to better characterize parton energy loss, we 
first build the px distributions for partons produced in 
given momentum bins. We then build the hadron px 
distributions that these partons originate and identify 
the corresponding momentum bins where these hadrons 
are produced. For the sake of definitcness, we carry out 
this procedure for exclusive momentum bins 2 — 4, 6 — 8, 
10 - 12 and 14 - 16 GeV. 

Figure Q] shows the parton and hadron pt distributions 
at central rapidity \r/\ < 0.5 for 2 — > 2 events, both in 
V + V ( e o = 0) and A + A (e = 2 GeV/fm) collisions at 
Y / J/vw = 200 GeV. On the left (right) column we show 
the leading (away) hadron px distribution. The leading 
hadron has been defined, event by event, as the hadron 
with the largest px, thus for each parton px bin, the 
away hadron momentum distribution is shifted towards 
lower values of px when compared to the leading one. 
The effect is more pronounced in A + A collisions due to 
energy loss. Notice also that momentum conservation at 
the parton level is reflected by the fact that the parton 
momentum distribution is the same for leading and away 
side hadrons. 

Figure [2] shows the parton and hadron px distributions 
for 2 — > 3 events, both in p + p (eo = 0) and in A + A 
(eo = 2 GeV/fm) collisions. In the figure, the left column 
shows the momentum distribution of the leading hadron 
which, as before is the hadron with the largest px in the 
event, the middle column displays the away side hadron 
selected as the hadron with the second largest momentum 
and the right column shows the away side hadron with 
the lowest px in the event. It is also interesting to note 
that the bulk of the hadron events with low momenta (in 
the range 0.2 GeV < px < 2 GeV), come from parton 
events with momenta in the range 2 GeV < px < 10 
GeV, whereas the hadron events with higher momenta 
(in the range 2 GeV < px < 4 GeV), come from parton 
events with momenta in the range 10 GeV < px < 16 
GeV. 

We now focus on the characterization of the events that 
contribute to the di-hadron azimuthal correlation either 
from 2 — > 2 or from 2—^3 parton processes. Figure [3] 
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FIG. 1: Parton (dotted lines) and hadron pr distributions for 2 — s> 2 events in p + p (en. = 0, solid lines) and in A + A 
(eo = 2 GeV/fm, dash-dotted lines) collisions for central rapidity \r)\ < 0.5 at ^/sjvjv = 200 GeV. On the 1 (2) column we show 
the pt distribution for the leading (away) hadron. The hadrons come from partons produced in four different momentum bins, 
from top to bottom: 2 - 4 (a), 6 - 8 (b), 10 - 12 (c) and 14 - 16 (d) GeV. 
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FIG. 2: Parton (dotted lines) and hadron pr distributions for 2 — s> 3 events in p + p (eo = 0, solid lines) and in A + A 
(eo = 2 GeV/fm, dash-dotted lines) collisions for central rapidity \r]\ < 0.5 at ^/sjvat = 200 GeV. From left to right we show the 
pr distributions for the leading (1) and the away hadrons (2, 3), which come from partons produced in four different momentum 
bins, from top to bottom: 2 - 4 (a), 6 - 8 (b), 10 - 12 (c) and 14 - 16 (d) GeV. 
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FIG. 3: Hadron pr distributions for events that contribute to the di-hadron azimuthal correlation. Shown are 2 — > 2 (a) and 
2 — > 3 (b) processes. The solid lines correspond to the case without medium and the dotted lines to the case within the medium. 
In both cases the leading hadron is in the momentum bin 2.5 GeV < pr < 4 GeV (1). For the upper row the away side hadron 
is in the momentum bin 1 GeV < pr < 2.5 GeV whereas in the second row, at least one of the away side hadrons is in the 
momentum bin 1 GeV < pr < 2.5 GeV (2). 



shows the hadron px distribution calculated under the 
aforementioned kincmatical conditions, where we use all 
the parton exclusive momentum bins, 2 GeV wide, in 
the range 2—18 GeV. The leading hadron is in the bin 
2.5 GeV < pr < 4 GeV. In the first row we show 2^2 
events where the away side hadron is in the momentum 
bin 1 GeV < pr < 2.5 GeV and in the second row we 
show 2 — > 3 events where at least one of the away side 
hadrons is in the momentum bin 1 GeV < pr < 2.5 GeV. 
We choose these hadron momentum bins to compare with 
the analysis of Ref. [l2j. Also, the selection of hadron 
events from parton pt bins in specified ranges favors 
mercedes-like configurations for hadrons. Notice that as 
expected, in each case the distributions without medium 
are consistently higher than the ones with medium ef- 
fects. 



V. THREE-HADRON AZIMUTHAL 
CORRELATIONS 

We now look at the azimuthal correlations for the away 
side considering the contributions from 2 — > 2 and 2 — s- 3 
hadron events. 

Figure HI shows the di-hadron azimuthal angular corre- 
lation normalized to the number of parton events and 
to the bin size from 2 — > 2 (upper row) and 2^3 
(lower row) events. For this figure we use parton px 
bins 2 - 4 and 10 - 12 GeV for p + p (e = 0) and A + A 
(e = 2 GeV/fm) collisions at ^/sWn = 200 GeV. Notice 
the typical peak shape at it radians in the 2 — > 2 case, 
that in our approach arises due to collincar fragmenta- 
tion, whereas in the 2 — > 3 case the peaks are located at 
27r/3 and 47r/3 radians and have a finite width. In the 
2 — 4 GeV parton pt bin, there is a marked difference 
between the height of the peaks for hadrons produced 
in the p + p and A + A cases. This is due to the fact 
that in the medium, the model considers a threshold at 
^thresh = 2.4 GeV below which the parton is unable to 



7 



1 

0.8 
2 0.6 
0.4 
0.2 



°o L 



(a1) 



— s„ = GeV/fm 
e„ = 2 GeV/fm 



2 4 
A(|> [rad] 



1 

0.8 
2 0.6 
0.4 
0.2 



if 



— e„ = GeV/fm 
e„ = 2 GeV/fm 



i [rad] 



0.05 
0.04 
j" 0.03 
| 0.02 
0.01 



(a2) 



. e„ = GeV/fm 
. e. = 2 GeV/fm 





> [rad] 



0.07 

0.06 

0.05 

1 0.04 

? 0.03 
z 

~ 0.02 
0.01 



(b2) e = GeV/fm 

e„ = 2 GeV/fm 

1 



2 4 
A(|> [rad] 



FIG. 4: Di-hadron azimuthal angular correlation normalized to the number of parton events and to the bin size from 2 — > 2 (1) 
and 2 — > 3 (2) events in parton pr bins 2 — 4 (a) and 10 — 12 GeV (b) for p + p (eo = 0) and A + A (eo = 2 GeV/fm) collisions 



hadronize and this value lies within the considered px 
bin. For all parton px bins higher than this one, the 
correlation of hadrons produced in both cases coincide, 
since the threshold value py resh is outside such bin. This 
is illustrated with the 10 — 12 GeV parton pr bin on the 
right column of Fig. 2) 

Experimental three-particle correlation functions are 
built from normalized distributions of a triplet of par- 
ticles chosen in a single event. Typically one of these 
three particles, the leading one, is chosen in a specified 
high transverse momentum bin while the two associated 
hadrons are chosen in a lower transverse momentum bin. 
By fixing the direction of the leading hadron, the anal- 
ysis is carried out in terms of the two angles that the 
away-side hadrons make with the leading one. Different 
correlation scenarios can lead to similar distributions. To 
distinguish between scenarios, it is customary to project 
the correlation choosing a certain combination of these 
angles to remain fixed and an other independent one to 
vary. 

One possibility is to look at the correlation function 
in terms of the angular difference A(f> = 0^ — 02 of the 
away side hadrons for a range of angles of one of the away 
side hadrons, say, 9% ■ Figure [5] shows this correlation in 
the proton and nuclear environments with ^/snn = 200 



GeV for a leading hadron momentum 2.5 GeV < pr < 4 
GeV and away side hadron momenta 1 GeV < pt < 2.5 
GeV, integrated over angular range 1.65 radians < 9Zf < 
2.2 radians. Shown are the normalized histograms for 
e = 2 GeV/fm (solid line) and e = (squares) ob- 
tained from the present Monte Carlo approach and from 
the approach in Refs. [|[ with eo = 2 GeV/fm (dotted 
line), compared to preliminary data from the PHENIX 
Collaboration [l2| . Notice that the correlation peaks 
at A(j> = 2.2 radians which is consistent with the in- 
terpretation that for the considered momentum bin, it 
is more likely that the two away-side hadrons are pro- 
duced close to a Mercedes-like configuration. This peak 
is also present in the data sample. Also, notice that the 
Monte Carlo approach histograms have a more defined 
peak than the approach where we regulated the diver- 
gences with angular cuts. The difference between the 
Monte Carlo result and the data away from A(j> = 2.2 ra- 
dians can be understood by recalling that we have imple- 
mented collinear fragmentation. To be able to populate 
the regions for A</> < 1 radians and Acf> > 2.5 radians, 
we would need to consider the contribution of the leading 
hadron and go beyond collinear fragmentation. 
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FIG. 5: Normalized di-hadron correlation in the nuclear en- 
vironment with eo = 2 GeV/fm as a function of the angular 
2 , integrated over angular range 1.65 



difference = 6»f 



2.5 GeV < pt < 4 GeV and away side hadron momenta 1 
GeV < pt < 2.5 GeV. The solid line and the squares corre- 
spond to the Monte Carlo approach with eo = 2 GeV/fm and 
eo = 0, respectively, and the dotted line to the approach in 
Refs. @| with e = 2 GeV/fm. Data are from the PHENIX 
Collaboration [121. 



VI. SUMMARY AND CONCLUSIONS 



fractions. Wc use the Monte Carlo approach to overcome 
the limitations set by collinear and soft divergences that 
otherwise have to be handled by means of angular cuts. 
We have studied the hadron px ranges that are populated 
when partons fragment. To be fully quantitative, one re- 
quires to generate Monte Carlo samples covering a wide 
range of parton px, both for 2 — > 2 and 2 — > 3 processes. 
Our approach on the other hand, based on populating 
hadron px bins from given parton px bins although re- 
ducing the need of a large Monte Carlo sample still has 
limitations, in particular, to assess quantitatively the rel- 
ative contributions between 2 — > 2 and 2 — > 3 processes. 
Work to overcome this limitation is underway and will 
be reported elsewhere. 

We conclude that since the energy loss model considers 
a threshold parton momentum py re,Bh = 2.4 GeV below 
which partons cannot fragment, there is a need to study 
the way those partons interact with the bulk medium and 
that do not hadronizc through fragmentation. These par- 
tons are rather abundant, since the low momentum part 
of their spectrum is more populated, and thus their effect 
must be considered. When these partons are produced 
in 2 — > 3 processes, due to momentum conservation they 
still come out from the hard scattering mainly with the 
characteristic mercedes-like configuration. Their interac- 
tion with the bulk partons should produce a transfer of 
momentum along these directions that will deform the 
medium producing showers of low momentum hadrons 
bearing the same distinctive angular distance of the orig- 
inal partons. This avenue is currently being pursued and 
the results will also be reported elsewhere. 



In this work wc have presented a Monte Carlo ap- 
proach to compute di-hadron correlation functions in- 
cluding contributions from 2 — > 2 and 2 — > 3 parton 
processes in p+p and A + A collisions. The production of 
parton events in p + p is implemented using MadGraph5 
where hadron events are obtained by evolving the partons 
into hadrons with collinear fragmentation by means of 
KKP fragmentation. For the nuclear environment we fol- 
low a modified fragmentation function approach whereby 
the hard scattering is assigned a random radial position r 
and a random direction <fi that the leading parton makes 
with the outward radial direction. These parameters are 
used to determine the average number of scatterings, the 
energy loss for each parton and their modified momentum 
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